arXiv:nucl-th/9811090vl 25 Nov 1998 


Lambda collective flow in heavy ion reactions 


Z.S. Wang, Amand Faessler, C. Fuchs and T. Waindzoch 

Institut fiir Theoretische Physik der Universitdt Tubingen, D-72076 Tubingen, 

Germany 


Abstract 

Collective flow of Lambda hyperons in heavy ion reactions at SIS energies is inves¬ 
tigated. It is found that a A mean field constructed on the basis of the quark model 
leads to a good description of the experimental data of the in-plane transverse flow 
of A’s. The attractive mean field can also give rise to an additional ’’virtual” A radial 
flow directed inwards, which is reflected by a ’’concave” structure of the transverse 
mass spectrum of the A hyperons emitted at midrapidity. The A radial flow is found 
to exhibit a strong mass dependence: The flow is visible in the Ni-|-Ni system, but 
is strongly reduced in the system of Au on Au. 
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1 Introduction 

Potentials of hyperons such as A, E, 5 etc. in a nuclear environment have 
attracted interest not only in nuclear physics, but also in astrophysics. It has 
long been realized that hyperons will appear in the dense core of neutron 
stars where the conversion of nucleons to hyperons through the weak interac¬ 
tion is energetically favored because at the Fermi surface the total energy of 
the nucleons will otherwise exceed the mass of hyperons[l,2]. The presence of 
hyperons signihcantly softens the equation of state of dense hadronic matter. 
This leads to the prediction of a lower neutron star mass, which is closer to as- 
trophysical measurements than previous predictions for hyperon-free neutron 
stars[2,3]. Other possible mechanisms to reduce neutron star masses include 
kaon condensation which was suggested by Kaplan and Nelson[4]. Hyperon- 
nucleus potentials can be extracted from hypernuclei[5]. Theoretical models 
have also been used to study hyperon-nucleus potentials. The approaches vary 
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from meson-exchange [6-8] to QCD (Quantum Chromodynamics)-based mod- 
els[9,10]. According to the quark model, a hyperon interacts with nucleons 
predominately through its non-strange quark content [11]. Consequently, po¬ 
tentials of different hyperons in nuclear matter should satisfy a simple scaling 
in terms of the number of non-strange quarks in the hyperons[10]. This spec¬ 
ulation should be confronted with precise experimental tests. 


Heavy ion reactions can serve as an alternative way for the investigation of 
hyperon-nucleus potentials. In fact, heavy ion reactions provide the only way 
to study in laboratories hyperon properties in a nuclear environment which is 
much denser and hotter than nuclei in their ground states. On the other hand, 
collective flow analysis turns out to be powerful for the study of the in-medium 
particle dynamics in heavy ion reactions. Experimental data about the flow 
of nucleons have been shown to yield stringent constraints on the nuclear 
equation of state (EOS) [12]. Recent studies revealed that the flow of kaons 
provides valuable information on kaon properties in the nuclear medium[13- 
15]. In this paper we will study the flow of Lambda hyperons in heavy ion 
reactions. 


We will explore in the present work two types of collective flows of A’s, namely 
the in-plane transverse flow and the radial flow. The transverse flow has been 
intensively studied for both nuclear fragments and for secondary particles, 
while the radial flow has been known conventionally as a collective expansion 
of the dense and hot hreball formed in colliding nuclei. In a recent work[15], 
we have shown that secondary particles such as and K~ mesons can also 
develop a collective motion in the same direction as the radial expansion of the 
fireball. This novel collective flow of the secondary particles is called ’’radial 
flow” as well. The present work will investigate the possible formation of the A 
radial flow. We assume in the present study heavy ion beam energies between 
1 and 2 GeV/nucleon. This energy region is near the A production threshold in 
free nucleon-nucleon collisions (1.58GeV), and is well studied experimentally 
at Berkerley and Darmstadt. This paper is organized as follows. In section 2, 
we describe the A dynamics in heavy ion reactions within the framework of 
Quantum Molecular Dynamics (QMD). Section 3 is devoted to the transverse 
flow of A’s. In section 4 we explore the radial flow of A’s. Section 5 is a 
summary of the paper. 


2 Lambda dynamics in the QMD model 


In heavy ion reactions at beam energies of 1-2 GeV/nucleon Lambda’s are pro¬ 
duced mainly via the reactions B-l-B —> B-l-A-l-iL’'' and tt-I-B ^ A-|-iL+[16]. 
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Here B stands for a nucleon or a nucleonic resonance A or N*. Due to 
strangeness conservation a A hyperon is accompanied by a positive charged 
kaon in both reactions. So we treat the Lambda production in the QMD model 
in a similar way as we did for mesons in a previous work[14]. We have used 
for both reactions calculated cross sections which agree with the experimental 
data[17,18]. While the baryon-baryon (B-B) reaction is mainly isotropic, the 
TT-baryon (vr-B) reaction exhibits a strong asymmetric angular distribution due 
to the P-wave contribution. The anisotropy of the ttB reaction has been taken 
into account with a parametrization of the experimental data. This enabled us 
to reproduce very well the experimental angular distribution of iP+’s in heavy 
ion reactions[19]. A Lambda hyperon can not be absorbed in nuclear matter. 
It, however, experiences strong elastic scatterings on nucleons. The AN elastic 
cross section is about 100 mb at a A laboratory momentum of 0.2 GeV/c, 
and becomes even larger with decreasing momenta[20]. We include Lambda- 
nucleon scatterings with use of a parametrization of the experimental cross 
section [21]. 


The A mean field in nuclear matter is constructed on the basis of the quark 
model as usual[ll,21]. The vector and scalar potentials acting on a A hyperon 
are taken to be two thirds of that experienced by a nucleon, namely 

= jE" (1) 

Ej- = jEj? (2) 

where and Ey are the nucleon scalar and vector potentials evaluated in 
the non-linear a-u model[22]. The mean-held potential of Lambda’s is then 
dehned as the difference between the in-medium dispersion relation and the 
free one: 

Uk = ^(rriA - E^)2 + p2 + s]) (3) 

In Fig.l we present the A potential at zero momentum as a function of nuclear 
matter densities. This potential agrees at saturation density po = 0.16/m“^ 
roughly with the value extracted from hypernuclei[23]. One can see from the 
hgure that the A potential is attractive at nuclear matter densities below 
three times the saturation density (p < 3po). The maximal density reached in 
heavy ion reactions at 1-2 GeV/nucleon is about 3po. Therefore, one will see in 
these reactions mainly effects of an attractive potential. We notice that the A 
potential (eq.3) is dehned in a non-covariant way since the contribution of the 
spatial components of the nucleon current is neglected[24]. This approximation 
should be justihed for the present study due to twofold reasons. First, the 
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velocity of A’s is quite limited in the center-of-mass frame of the colliding 
nuclei, since we consider in the present paper a beam energy close to the A 
production threshold for free NN collisions. Second, the A hyperons experience 
frequent scatterings with nucleons, which further reduce dramatically the A 
velocity relative to surrounding nucleons. Therefore, the flow of the Lambdas 
should not change by much as one changes from the non-covariant A dynamics 
to a fully covariant one. 


The motion of A hyperons is determined with the Hamiltonian: 

+ (7,A . (4) 

*A ^ '' 


The A production probability in an in-medium hadron-hadron collision could 
be different from that of a collision in free space. The magnitude of this 
medium effect depends on the in-medium potentials of all the hadrons in¬ 
volved in the collision, i.e., nucleons, pions, kaons and A's. The potentials of 
these hadrons have been not yet well determined both experimentally and the¬ 
oretically. Uncertainties in evaluating the medium effect may also arise from 
the procedure to decompose the interactions of the hadrons with the nuclear 
medium into a long-range mean-field part and a short-range hard-scattering 
part. The magnitude of the medium modification of the A and produc¬ 
tion probability is still under current debates[17,25]. However, such an effect 
should not change the collective flow of the A’s, and therefore is neglected in 
the present work. 


3 Lambda transverse flow 


Transverse flow is a collective emission of particles sidewards in the reaction 
plane of colliding nnclei. This type of collective flow can be qnantitatively 
described with the particle momentnm projected onto the reaction plane as 
a fnnction of rapidity[26]. A flow results in a non-zero in-plane momentnm 
at projectile and target rapidity. The existence of the transverse flow was 
originally predicted for nnclear matter by the hydrodynamic model[27,28], and 
was later verifled by experiments[29]. It has been fonnd that the transverse 
flow of nncleons or composite particles snch as denterons, helinms etc. exhibits 
interesting dependence on the bombarding energy. At low energies where the 
attractive mean held dominates the interaction, the projectile fragments are 
deflected towards the target [30]. If the bombarding energy is high so that the 
overlapping matter is strongly compressed, positive pressure will develop in 
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the interaction zone. Consequently, the projectile fragments will be deflected 
away from the target (so-called ”bounce-ofF’)[29]. Between the two extreme 
cases, one hnds a certain incident energy where the attractive component and 
repulsive component of the nuclear mean held counterbalance each other, and 
therefore the transverse how disappears[31-33]. This special energy, called in 
literature ’’balance energy”, is about Ebai = 10-25 MeV/nucleon as found 
experimentally[31,32]. The beam energies considered in the present work (1-2 
GeV/nucleon) are much higher than E^ai- We therefore will see a bounce-oh 
of nucleons for the reactions under consideration[14]. 


The studies of the transverse how of nuclear fragments have motivated a lot 
of recent studies which analyzed the in-plane momentum of secondary parti¬ 
cles such as pions, kaons, antiprotons etc. at varying rapidity[13,14,34,35]. In 
analogue to the nucleon how, one usually says that a transverse how of the sec¬ 
ondary particles happens as a non-zero in-plane momentum of the particles 
is observed at projectile or target rapidity. The transverse how of the sec¬ 
ondary particles turns out to provide information on the in-medium dynamics 
of these particles. The pion transverse how was found to be anticorrelated to 
the nucleon how[34,35]. This behavior of the pion how can be attributed to 
rescatterings and reabsorptions of the pions in the nuclear medium. The kaon 
transverse how was shown to be sensitive to in-medium modihcations of kaon 
properties[13,14]. In this section we explore the transverse how of A’s, and its 
possible sensitivity to the Lambda-nucleus potential. 


We present in Fig.2 the A in-plane momentum as a function of rapidity cal¬ 
culated from the QMD model for the reaction Ni-|-Ni at an incident energy of 
1.93 GeV/nucleon and within an impact parameter region of b = 0-4 fm. Also 
shown in the hgure are the experimental data of the FOPI Collaboration[36]. 
A transverse momentum cut of Pt/m\ >0.5 has been used in obtaining both 
the experimental and the theoretical spectrum. A rapidity is calculated in the 
center-of-mass frame of the reaction system and normalized to the projectile 
rapidity. We have performed the QMD calculations in two cases: with and 
without the in-medium A potential. One can see from the hgure that the cal¬ 
culation without the potential deviates from the experimental data at high 
rapidity, while the full calculation reproduces well the data. The FOPI ex¬ 
periment has found that the A how essentially follows the proton how in the 
same reaction, while the positive charged kaon how is very weak (The kaon 
in-plane momenta are nearly zero.). In a previous study with the same QMD 
model as used in the current work, we have got a very good agreement with 
the FOPI data for both proton and kaon how[14]. As can be seen in Fig.2 
the A how obtained without the potential is weaker than the one from the 
full calculation. One can understand the ehect of the potential on the A how 
in a similar way as one did for the kaon how. Provided that the Lambda’s 
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and kaons would experience neither rescattering nor an in-medium potential, 
they should have a same flow pattern, since a A hyperon and a meson are 
created simultaneously from a BB or ttB collision. The mesons feel in the 
nuclear medium a repulsive potential which repels the kaons away from their 
sources and by that leads to a near-vanishing transverse flow. A hyper- 
ons, on the contrary, experience an attractive mean held at the beam energies 
under consideration. Thus the A flow changes to be more pronounced due to 
the potential. 


Aside from the A mean held, A-nucleon scatterings also affect the A flow. 
At the beam energy under consideration the projectile-like and the target- 
like remnants are deflected away from each other as mentioned above. The 
transverse motion of these remnants can transfer to the Lambda’s through 
the A-nucleon scatterings. Consequently, the A flow is enhanced due to the 
A-nucleon scatterings. 


The A transverse flow of the same reaction has also been studied in Ref. [21] 
with a different transport model, i.e. the Relativistic Boltzmann-Uehling- 
Uhlenbeck (RBUU) model. Authors of Ref. [21] used a non-covariant descrip¬ 
tion of the Lambda dynamics which is quite similar to the one adopted by 
the present paper. Good agreement is found between Ref. [21] and the present 
study concerning the transverse flow of iL+’s and A’s. As far as the proton 
flow is concerned, small discrepancies exist between the two theoretical stud¬ 
ies. Our QMD model turns out to reproduce very well the proton flow data up 
to ycm/Uproj = 1-5, while the RBUU model of Ref.[21] slightly underpredicts 
the data at ycm/yproj ^1* 


4 Lambda radial flow 


Radial flow conventionally means a subsequent expansion of the dense and 
hot fireball formed in heavy ion reactions[37,38]. The fireball expansion has 
similar origin as the bounce-off: both are driven by the positive pressure in the 
overlap region of colliding nuclei. The threshold beam energy for onset of the 
expansion is found experimentally to be about 35 GeV/nucleon[39], slightly 
higher than the threshold energy for the occurrence of the bounce-off phe¬ 
nomenon, Ehai- The fireball expansion gives rise to some non-thermal behavior 
of nucleons[37], composite particles such as deuterons, tritons etc.[40]. It can 
also implement non-thermal features to secondary particles, provided that the 
beam energies are so high that hadron-hadron collisions in the late expansion 
phase of the reaction are still energetic enough to produce these particles. 
This condition is fulfilled at AGS and SPS energies ( about 10 GeV/nucleon 
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and 150 GeV/nucleon, respectively) for pions, kaons, Lambda’s etc.. There¬ 
fore, one can find information on the fireball expansion not only from nucleon 
observables but also from observables of these secondary particles[41-44]. At 
a subthreshold beam energy, one will hardly hnd effects of the hreball ex¬ 
pansion on secondary particles, since the particles are produced mainly in 
the compression stage prior to the subsequent hreball expansion. In a recent 
study, we demonstrated that and K~ mesons produced at a subthresh¬ 
old beam energy (1.58 GeV/nucleon for iL+’s and 2.5 GeV/nucleon for iL“’s) 
can exhibit non-thermal features, however, due to their in-medium potentials 
instead of the hreball expansion[15]. ( K~ ) mesons feel a repulsive (at¬ 

tractive) potential in nuclear matter, which mutually accelerates (decelerates) 
the {K~) mesons as they escape from the dense hreball. This results in 
a ”shoulder-arm” (’’concave”) structure in the transverse mass spectrum of 
the {K~) mesons emitted at midrapidity, which obviously deviates from a 
Boltzmann distribution. We called kaon radial how the common acceleration 
or deceleration by the in-medium potentials. 


In this section we investigate if A hyperons can also develop a radial how 
similar to kaons. Fig.3 shows the transverse mass spectrum of the A hyperons 
emitted at midrapidity ( -0.4:<ycm/yproj< 0.4 ) calculated for the same reac¬ 
tion as in Fig.2. A ’’concave” structure can be clearly observed in the spectrum 
at small values of the transverse mass (m^- Ma = m\ + p\-m\ < 0.2 GeV). 
This feature obviously distinguishes the A spectrum from a Boltzmann distri¬ 
bution (see a Boltzmann ht to the high-energy part of the A spectrum which 
is also shown in Fig.3). The ’’concave” structure of the A spectrum arises from 
the attractive A mean held. If the A potential is neglected, the ” concave” shape 
disappears. This can be seen in Fig.4 where the A transverse mass spectrum 
from the QMD calculation without the potential is presented. The A spectrum 
is now exactly a Boltzmann distribution. 


In order to extract information on the in-medium potentials from the trans¬ 
verse mass spectrum, we suggested in Ref. [15] an equation to ht the spectrum. 
The ht function has been derived by boosting a thermal Boltzmann distribu¬ 
tion with a common radial velocity due to the in-medium potential. 


d^N 


+ ^aT{^ + 1 ) + 




d(j)dymtdmt 


^{jE + aTy 




pz. 




P 


-( 5 ) 


where E = nitcoshy, p = y/p^ + mysinh^y, a = 'yjdp/T, 7 = (1 — 

There are two free parameters appearing in the ht equation (eq.5). Both of 
them have clear physical meanings: T is the average temperature of A sources, 
while the parameter (3 = v/c is a common radial velocity of the Lambda’s 
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caused by the potential. An attractive potential leads to a reduction of A 
radial velocities. We thus should use in the £t equation a f3 directed towards 
the center of the fireball. One can extract the T value from the A spectrum 
which comes from the QMD calculation without the A potential. This yields 
a value of T = 115 MeV. The (3 value can be evaluated by comparing the 
simplified QMD calculation with the full one including the potential. We found 
from the calculations that the average mass of the midrapidity Lambda’s is 
reduced by the attractive potential from < rriT >-mA = 127 MeV to < uit >- 
m\ = 115 MeV. This reduction accounts for a f3 value oi {3 = 0.025 in the 
fit procedure. In Fig.3 we plot the result of the fit by using eq.5 with T and 
(3 parameters determined as above. A satisfactory agreement can be found 
from the figure between the full QMD calculation and the fit procedure. This 
demonstrates that a A radial flow forms in the reaction due to the attractive 
potential. However, the energy reduction caused by the potential is rather 
small compared to the thermal energy: the former is only about 10% of the 
latter. Therefore, experimental data of high accuracy are required in order 
to identify the A radial flow. We have noticed that the EOS Collaboration 
reported a measurement of the transverse mass spectrum of midrapidity A’s 
in the reaction Ni + Cu at 2 GeV/nucleon[45]. However, the error bars of the 
EOS data are too large (in particular in the region of 0.2 GeV< my-mA < 0.5 
GeV) to enable us to draw any conclusion about the A radial flow. 


It is interesting to study the influence of the size of the reaction system on 
the A radial flow. In order to do this, we also study with the QMD model 
reactions induced by massive nuclei, i.e. Au+Au. Fig.5 shows the calculated 
transverse mass spectrum of midrapidity Lambda’s for this reaction at an 
incident energy of 1 GeV/nucleon and at an impact parameter of b = 3 fm. 
This spectrum can be described very well by a pure Boltzmann distribution 
as can be seen in the figure. In Fig.6 we present the A spectrum of the same 
reaction as in Fig.5 from a calculation where the A mean held is neglected. 
One hnds again that the QMD result without the A potential agrees well with 
a thermal distribution. The A spectrum from the simplified calculation is very 
similar to the one from the full calculation (Fig.5). The average transverse 
mass < rrit > is only reduced by about 1 MeV due to the potential. Thus, no 
effect of the A mean field can be observed and the radial flow is invisible in 
the Au+Au reaction. 


That A radial flow forms in the Ni+Ni reaction but disappears in the Au+Au 
reaction can be attributed to A-nucleon scatterings. It is clear from Fig.4 and 
Fig.6 that A-nucleon scatterings alone can not induce any two-temperature 
structure of the A spectrum. However, as soon as a A radial flow develops as a 
consequence of the attractive potential, A-nucleon scatterings begin to play a 
role in canceling the effect of the potential, since the scatterings bring the A’s 



from the low energy region of rriT - rn\ < 0.2 GeV to higher energies. As we 
have mentioned, a A hyperon has a short mean free path in nuclear matter, 
especially for low-energy Lambda’s. One expects more A-nucleon scatterings 
in a larger system, and thus less low-mr Lambda’s. Consequently, one hnds 
no ’’concave” structure in the A spectrum in the reaction Au-|-Au, while this 
structure survives in the Ni-|-Ni reaction as a result of reduced A-nucleon 
scatterings. The density dependence of the A potential could also play a role 
in the A radial flow. As can be seen in Fig.l, the A potential changes to be 
shallower as the density increases beyond the normal nuclear matter density 
Po- A larger system has a larger stopping power, and thus a longer lifetime of 
the dense hreball. Consequently, the A hyperons experience a less attractive 
potential in a larger system. This also leads to a weaker radial flow in the 
Au-|-Au reaction than in Ni-|-Ni reaction. 


It is worthwhile to indicate that the A radial flow is very sensitive to inci¬ 
dent energies. The conclusions of the present work are drawn for the incident 
energies close to the A production threshold for free NN collisions. These rel¬ 
atively low beam energies, i.e. 1-2 GeV/nucleon, lead to a negligible effect 
of the fireball expansion on the A’s as we have mentioned. However, as one 
goes to AGS energies or even SPS energies, the hreball expansion will play 
a non trivial role in determining the A spectrum. The hreball expansion can 
reduce the multiplicity of low-mr A’s by pushing them to a hnite niT corre¬ 
sponding to the common expansion velocity of the hreball. Consequently, the 
transverse mass spectrum of the A’s will show a ” shoulder-arm” shape rather 
than a ’’concave” shape. Such ”Shoulder-arm” structure has been observed by 
the E891 Collaboration for the incident energy of 11.6 GeV/nucleon[43]. This 
observation does not contradict with the present study. 


5 Summary 


In this paper we have studied the collective how of A hyperons in heavy ion 
reactions at SIS energies with the QMD model. It is shown that the FOPI 
data of the A transverse how, which is quantitatively described with the A 
in-plane momentum as a function of rapidity, can be well reproduced by using 
a A potential constructed based on the quark model. The QMD model used in 
the present paper can also describe very well the measured proton and kaon 
how as shown in our previous work. A radial how of A hyperons is found to 
arise from the attractive A potential in the Ni-|-Ni reaction. The how results 
in a ’’concave” structure of the transverse mass spectrum of the A hyperons 
emitted at midrapidity. However, as one changes to the reaction Au-|-Au, the A 
radial how vanishes. The size dependence of the A radial how can be attributed 
to frequent A-nucleon scatterings and reduced attraction of the A potential in 
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a large system. 
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Figure Captions 


Fig. 1. The A mean-field potential in nuclear matter used in this work (at 
zero momentum). The circle denotes the potential at normal nuclear matter 
density po extracted from hypernuclei experiments[23]. 


Fig. 2. The in-plane A momentum as a function of rapidity from the Ni-|-Ni 
reaction at 1.93 GeV/nucleon and within an impact parameter region of b=0-4 
fm. The solid line denotes the QMD calculation including the full in-medium 
Lambda dynamics, while the dashed line is the QMD calculation without the 
A potential. The experimental data from the FOPI Collaboration are also 
presented. 


Fig. 3. The transverse mass spectrum of the A hyperons emitted at midrapidity 
( -0.4<?/c.m./l/proj<0.4 ) from the same reaction as in Fig.2. The histogram 
presents the result of the QMD calculation which is performed including the 
full in-medium Lambda dynamics. The dashed line presents the result of a 
pure thermal Boltzmann fit to the high energy part of the spectrum ( rrit-mo 
> 0.2 GeV ), while the dotted line is the fit according to eq.(5) where one 
assumes a A radial flow in addition to the thermal motion. 


Fig. 4. The transverse mass spectrum of the A hyperons emitted at midrapidity 
( -0.4<|/c.m./|/proj<0.4 ) from the same reactions as in Fig. 2. The histogram 
presents the result of the QMD calculation without the A mean field. Now a 
pure thermal fit (lines ) is sufficient to reproduce the spectrum. 


Fig. 5. The transverse mass spectrum of the A hyperons emitted at midrapidity 
( -0.4<|/c.m./|/proj<0.4 ) from the reaction Au-|- Au at 1 GeV/nucleon and b=3 
fm. The histogram presents the result of the QMD calculation with the A mean 
field, while the line is a Boltzmann fit to the QMD result. 


Fig. 6. The transverse mass spectrum of the A hyperons emitted at midrapidity 
( -0A<yc.m./yproj<0.4: ) from the same reaction as in Fig.5. The histogram 
presents the result of the QMD calculation without the A mean field, while 
the line is a Boltzmann fit to the QMD result. 
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